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A general theory of molecular-scale energy transfer promises to yield insights into the mechanisms 
of molecular machines and to unify a diverse range of experiments and phenomenological models. 
We progress towards a general theory by deriving analytic solutions to Brownian motion on a 
multi-dimensional tilted periodic free-energy potential. Analogous to the tight-binding model of 
quantum mechanics, we consider the limit of deep potential wells and show that the continuous 
theory transforms to a discrete master equation describing hopping between localized (meta-)stable 
states. For non-separable potentials the master equation describes energy transfer between degrees 
of freedom and we derive expressions for the efficiency and rate of energy transfer valid beyond the 
near-equilibrium limit. Our predictions are consistent with non-equilibrium fluctuation theorems 
and the Onsager relations and provide an opportunity to experimentally test the theory. 
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Biological cells use specialized proteins to convert and 
utilize chemical energy while operating far from equilib- 
rium, with minimal inertia, and in the presence of signif- 
icant thermal fluctuations [TH3]- Insights into the mech- 
nisms of these molecular machines are being provided 
by single-molecule experiments [TJ [2 and the arti- 
ficial synthesis of molecules that mimic motor proteins 
[3] 15HTT] . Building on one-dimensional theoretical stud- 
ies of stochastic non-equilibrium transport 2 , Brownian 
motion on a multi-dimensional free-energy landscape has 
been proposed as a mechanism for molecular-scale energy 
transfer [HI [13]. 

In this paper we begin with the continuous stochas- 
tic theory for overdamped Brownian motion on a multi- 
dimensional tilted periodic potential and derive a dis- 
crete master equation valid in the limit of deep poten- 
tial wells. This master equation describes hopping be- 
tween localized (meta-)stable states and can be solved 
analytically. We demonstrate that this approach has the 
key elements required for a general theory of molecular- 
scale energy transfer: (i) it is based on a physically com- 
pelling microscopic theory, (ii) it is consistent with well- 
established non-equilibrium thermodynamic results such 
as non-equilibrium fluctuation theorems |14[ I15j and the 
Onsager relations [TBI E]> an d (hi) it provides specific 
testable predictions. 

We consider overdamped Brownian motion on a multi- 
dimensional tilted periodic potential described by the 
Smoluchowski equation |18j 

^M =£/P(M) = _ V . J(r,t), (1) 

where P(r, t) is the probability density of finding the sys- 
tem at position r at time t and J(r,t) is the probability 
current density. Each degree of freedom is a generalized 
coordinate capturing the main conformal motions of the 
molecules and representing displacements in real space 
or along reaction coordinates [15] . In Eq. 0, the evolu- 
tion operator is defined by jCf = k B TV 2 + V • VV^(r), 



where T is the temperature, ks is Boltmann's constant, 
and 7 is the damping constant. For notational simplicity 
we choose uniform damping although the generalization 
to coordinate variations is straight forward |20j . 

The free-energy potential Vf(r) has both entropic and 
mechanical contributions [13]. We assume a potential in 
the form of a rapidly varying periodic part and a slowly 
varying linear tilt [13], i.e., Vf(r) = Vo(r) — f ■ r, where 
V (r+An) = V (r), Ajj> — djSjji and n = J2j n j^j with 
rij integer. The minima of the periodic potential (cen- 
tered at r = An) define microscopic (meta-)stable states 
and the potential barriers between minima determine the 
probability of hopping transitions. The linear potential 
represents a constant macroscopic thermodynamic force 
due to an external mechanical force or an entropic force 
such as a concentration gradient across a membrane or 
an out-of-equilibrium chemical concentration. 

In the / = case, the steady state has the form 
P ss (r) cx exp[—Vo(r) / ksT] and the probability current 
density vanishes everywhere, i.e., J S s(f) — 0. In con- 
trast, for f j^z the force induces a finite current and 
drives the system out of thermal equilibrium [2] . Energy 
transfer occurs when the force in one coordinate gives 
rise to current in another. This is only possible when the 
potential Vf(r) contains a non-separable term [12] . 

The above formalism provides a from-first-principlcs 
mathematical framework that encompasses all molecular- 
scale energy transfer, including energy conversion in cy- 
toskeletal motors, rotary motors such as ATP synthase, 
and ion pumps [4]. This theory also provides a com- 
pelling physical picture of a molecular-scale system un- 
dergoing Brownian motion on a potential landscape that 
directs the average behavior of the system thereby en- 
abling energy coupling between degrees of freedom. Even 
in the two-dimensional case, the Smoluchowski equation 
([I]) is not analytically tractable in general and is difficult 
to connect with experiments, phenomenological models, 
and other general results from non-equilibrium thermo- 
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dynamics. In the experimentally relevant regime of deep 
potential wells, we transform Eq. ([I]) to a master equa- 
tion that can be solved analytically. 

Analogous to the tight-binging model of a quantum 
particle in a periodic potential [211 , we consider that 
the periodic potential Vo(r) creates deep potential wells. 
The microscopic (meta-)stable states centered at each 
potential minima are strongly localized and the evolu- 
tion operator Cf induces infrequent hopping transitions. 
Formally, VV/(r) is periodic so we invoke Bloch's the- 
orem [21 . In the / = case, the eigenvalue problem 
for EqTjj} has the form C o 0° k (r) = -A° fc < fc (r), 

where $£fc(r) = e* fc-r u° fc (r) and u° k (r) has the pe- 
riodicity of the periodic potential. As for a quantum 
particle in a periodic potential, the eigenvalues A° k are 
real and separate into bands denoted by the band in- 
dex a and k is confined within the first Brillouin zone 
[21]. Periodic boundary conditions are used with pe- 
riod Na where AT is a diagonal matrix of integers. 
We choose an infinite spatial extent (Njj — > oo) which 
means k is continuous. Similar results can be obtained 
for rotary motors by taking Njj finite and k quan- 
tized. The eigenstates 0° k (r) form a complete orthonor- 
malset |2Q]: ^(rlMEJg^W^M = 
5(r-r>) and J dr exp[U o (r)/fc s T]0°* fe (r)0°, fc ,(r) = 
5 aa '8(k — k'), where the B subscript on the k integral 
denotes a single Brillouin zone. Unlike for a quantum 
particle, A° k > is to be interpreted as a decay rate 
rather than a frequency/energy and the ground state is 
#} j0 (r) oc exp[-F (r)/fc B T] with X° = 0. 

In the tight-binding regime it is convenient to 
define the localized Wannier states w a . n (r) = 
D f B dk<f>° k (r)ex.p(-ik- An), where D = Y[j(aj/2ir). 
These Wannier states also form a complete orthonormal 
set and we expand the probability density as P(r, t) — 
(1/ D)J2 a ,nPa,n{t)w atn (r). Using this expansion, Eq. 
([!]) can be transformed to a discrete form. In general 
this discrete master equation can be interpreted in terms 
of hopping between Wannier states within a band and 
coupling between bands induced by the tilting potential. 
Assuming that the tilt is sufficiently weak that inter-band 
coupling is small compared to gaps in the band struc- 
ture, inter-band coupling can be neglected. In this case 
the band structure of eigenvalues enables a separation 
of timescales between the rapidly decaying higher bands 
and the slowly evolving lowest band governing the long- 
time behavior of the system. Retaining only the lowest 
band and dropping the band subscript for the remainder 
of this paper, we write the resulting master equation [20] 



dp n {t) 
dt 



EK-nA'W-K^JnW], (2) 



where 



k£ = dre v °^l kBT w n {r)CfW {r). (3) 



In the tight-binding limit where the system is stongly 
localized at the minima of the potential wells, the sum- 
mation in Eq. ([2| need only be extended over nearest 
neighbors. The coefficients p n (t) are real and can be in- 
terpreted as the probability that the system is localized 
in the nth potential well. The hopping rates k£ are also 
real and satisfy J2n K n = 0- If the potential Vo(r) is 
separable, the hopping rate ^ has the form 



^2 k ™ 



ra,/0, 



(4) 



and hopping transitions occur independently along each 
coordinate. Energy transfer can not occur in this case. 

The hopping rate k£ is in general a complicated func- 
tion of /. However, in the case of tight binding and 
weak tilt a simple functional dependence can be derived. 
The Smoluchowski operator can be written in the form 
C f = e^p[-Vf(r)/2k B T]H f exp[Vf(r)/2k B T], where 
7%/ = k B TW 2 — Uf(r) is a self-adjoint operator with 
the potential 17/ (r) = [VV (r) - f} 2 /Ak B T -V 2 V (r)/2 
18 . In the tight-binding approximation when the tilt is a 
small perturbation near the potential minima, the Wan- 
nier functions decay fast compared to the slowly varying 
linear potential and we write w n (r) exp(/ • r/2k B T) » 
w n (r) exp(/ • An/2k B T). Furthermore, the main contri- 
butions to the integral ([3| come from regions around the 
minima of Uf(r) where the localized Wannier states are 
non- vanishing. Assuming the tilt is sufficiently small that 
it does not significantly alter the position of the / = 
potential minima, we take Uf(r) w Uo(r). Combining 
these results yields [2U] 



1 

D 



dr e Vo{r)/kBT w n {r)e 



-V (r)/2k B T^ V (r)/2k B T 



f-An/2k B T 



xe 

f-An/2k B T 



w (r) 



where 



D I dk X° k e lk An 



(5) 



(6) 



We estimate the validity regime for Eqs. Q and ^ 
by approximating the potential in the vicinity of its min- 
ima by a harmonic potential. In this limit the Wannier 
states are the harmonic oscillator states of the potential 
minima and the separation between the lowest and first 
eigenvalue bands is of order fl, the smallest eigenvalue 
of the matrix Mjj> = d 2 Vo{r)/drjdrji\ r —An/l- There- 
fore, all bands except the lowest can be neglected for 
t > 1/fi. To derive Eq. we have also assumed 

that f 2 /2k B T -C 7f2. Under this condition the tilt is 
small enough that transitions to higher eigenvalue bands 
of the / = potential are negligible. When a single 
maximum of Vb(r) separates the potential minima, the 
hopping rates k£ can be approximated by Kramer's rela- 
tion 20, 22, 23J. In fact, the physical justification for Eq. 



3 



^ closely parallels the physical argument in the deriva- 
tion of Kramer's relation (231 124] : rapid relaxation within 
potential wells accompanied by slow transitions between 
wells. The hopping rates k° in general can be calculated 
directly from the free-energy potential of the system and 
the potential can be determined by single-molecule ex- 
periments |25j or molecular dynamics simulations (e.g., 
[26]). 

The transformation from a continuous stochastic 
multi-dimensional formalism to a discrete master equa- 
tion represents a significant simplification of the system 
dynamics and has been attemped elsewhere [Jj2 l?7H2"9"] . 
The benefit of our approach is that we make this trans- 
formation explicit providing analytic expressions for the 
hopping rates. Furthermore, in the tight-binding limit, 
Eq. f|) gives the functional dependence of the hopping 
rates on the macroscopic thermodynamic force. Solving 
the master equation ^ to determine physical properties 
of the system provides an opportunity to test the the- 
ory against established non-equilibrium thermodynamics 
results and experiment. 

General non-equilibrium thermodynamics results are 
not imposed as constraints on our theory but instead can 
be determined analytically. In particular, the hopping 
rates ^ can be used to determine the ratio of forward 
to backward hopping, i.e., 

c/ 

(7) 



_ e f An/k E 



Equation ^ is consistent with general non-equilibrium 
fluctuation theorems [HI [T3] . 

The spatial drift v = d(n)/dt — J2 n ndp n (t)/dt repre- 
sents the average rate of hopping and can be determined 
from the master equation to be 



(8) 



Equation ([8]) shows the functional dependence of the drift 
on the macroscopic thermodynamic force, and vanishes 
for / = 0. Interpreting Xj = fjCij/T as the gener- 
alized thermodynamic forces and Vj as the conjugate 
fluxes, Eq. ^ represents a generalized force-flux rela- 
tion. Near equilibrium, \Xj\/kB "C 1 and Eq. ^ re- 
duces to v — J2n nK nJ2j XjUj/2kB- Therefore, vj — 



which 



Ej' L J3' X J' with L n' = 'E n n 3 n j' K r l / 2k B 

satisfies the Onsager relations [TBI E] ■ 

Force and velocity are directly measurable in single 
molecule experiments [7J H31 [30] . To explore these prop- 
erties in detail, consider the conceptually simpler two- 
dimensional case. The force-flux relation (JsJ) becomes 

V x = ( K fi. ) _ K f-1,0)) + ( K fl,l) _ K (- 



where in the nearest-neighbor sum we have assumed that 



and 



"(0,1) 



(i,o)U K (o,i)l'l K (i,i)l- The no PP in g rates 
represent transitions occuring indepen- 



dently along each coordinate while the hopping rates 
K ?i l) represent coupling transitions that transfer energy 
between the two degrees of freedom. 

In the case that the linear potential is downhill in di- 
rection x (X x > 0) and uphill in direction y (X y < 0), 
energy transfer from x to y is thermodynamically viable 
when the coupling transitions are downhill (X x + X y > 
0). The efficiency of energy transfer can be determined by 
the ratio of the power output P ou t = —v v X y T to input 
Pin — v x X x T , i.e., 



Pn 



V 



p. 

1 in 



v y X v 
v x X x 



(10) 



Equation ( 10 1 can be written explicitly in terms of Xj by 
inserting Eqs. ([9|. The efficiency satisfies < r\ < 1, and 
Tj — > 1 when the independent hopping transitions are neg- 
ligible, |K° 0)1) |, |«( 1)0 )l < l«(l,l)l) and when x x -> -Xy 
The latter condition can be recognized as the approach to 
thermal equilibrium along the coupled coordinate. Equa- 
tion ( 10 1 can be interpreted as a trade-off between the 
power output P out and the efficiency r\ [30 . This trade off 
may have important biological consequences |31l 132] and 
its characteristics beyond the near-equilibrium regime 
will be explored in future work. 

The eigenvalue band structure plays a key role in 
determining the system properties. The master equa- 
tion ([2]) can be transformed to the diagonal form 
dck(t)/dt = — A^Cfe(t), where the eigenstates are Ck(t) = 
J2nP™(t) e ~ lk An an d the eigenvalues are 



= 4k° 10) sin (k x a x /2) sin (k x a x /2 -\ 
+4k° ^ sin (k y a y /2) sin (k y a y /2 
+4«( 1)1 j sin (k x a x /2 + k y a y /2) 
x sin (k x a x /2 + k y a y /2 + iX x /2ks 



iX x /2k B ) (11) 
f iX y /2k B ) 



lXy/2k B ) . 



-1,-11 



(9a) 



Equation (111 defines the lowest Bloch band in the 
tight-binding limit. The gradient of the imaginary part 
at the origin is proportional to the drift, i.e., v on 
VfcIm(Aj[)|fc = o, and the curvature of the real part at 
the origin is proportional to the time derivative of the 
covariance matrix [33 , i.e., dUjiiUj) — (rii) (rij)) / dt oc 
d 2 Re(Xl.)/dkidkj\k = o. Figure 111 shows contour plots of 
the real and imaginary parts of the eigenvalues through- 
out the first Brillouin zone for (a) weak coupling near 
equilibrium, (b) strong coupling near equilibrium, and 
(c) strong coupling far from equilibrium. From (a) to 



(b), the drift goes from v y < to v y > despite the fact 
2K& n smh(X x /2ks) + 2k% ^ sinh(X x /2fcs + X y /2fcs)the X y < 0. Only quantitative differences are observed 
J \j_r„f _J \ ( 9b l ar from equilibrium. 

To provide a concrete two-dimensional example, — 



( K (o,i) K fo,-i)) + ( K fi,i) 



2k° (oa) smh(X y /2k B ) 



2k °, 



( — 1; — l)' 

\ smh(X x /2kB + X y /2kB)$ider a coupled chemical reaction system composed of the 
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-2 



FIG. 1. (color online). Contour plots of the (upper) real and 
(lower) imaginary parts of \f k fe n/K( ^ (arbitrary units). 
The symbols + and — denote maxima and minima, respec- 
tively. The plots (a r ) and {at) correspond to X x /ks = 0.1, 
k° m) = 0.01k° 0i1) ; (b r ) and (6<) to X x /k B = 0.1, K° 1A) = 



5 k 



{01) ; and (c r ) and {a) to X x /ki 



Other parameters are Xy/ks = —0.05 and (tn 



v (i,i) 



5 k 



(o,i)- 



(1,0) 



"(o,i)- 



three elementary reactions 



A^B, C^D, A + C^B + D, 



(12) 



numbered 1 to 3 from left to right. The Gibbs free en- 
ergies AGj are the thermodynamic forces driving the 
system and thermodynamic consistency requires AG\ + 
AG2 = AG3. The net rate for each chemical reaction is 



R b i=R f Al 



,AG 3 /fc B T 



)• 



(13) 



where and are the forward and backward reaction 
rates, respectively, given by the usual mass-action expres- 
sions in terms of species activities and reaction rate con- 
stants. In our formalism, the generalized thermodynamic 
forces are Xj — —AGj/T and the generalized fluxes are 
Vx = r% + r3 and v y — T2 + r§. If reaction 1 is spon- 
taneous (AGi < 0), and reaction 2 is non-spontaneous 
{AG 2 > 0), reaction 3 enables energy transfer between 
reactions 1 and 2 and this occurs spontaneously when 
AG3 < 0. Reactions 1 and 2 represent dissipative leak 
reactions that by-pass the coupling mechanism of reac- 
tion 3 [33]. 

The benefit of a general theory of molecular-scale en- 
ergy transfer is the opportunity to gain deeper theoret- 
ical insights into the mechanisms and characteristics of 
molecular machines and to unify a diverse range of exper- 
imental results and phenomenological models. Brownian 
motion on a multi-dimensional tilted periodic potential 
provides a compelling candidate for such a general theory. 
On that basis, we have derived a discrete master equa- 
tion in the limit of deep potential wells. This master 
equation is consistent with non-equilibrium fluctuation 
theorems and the Onsager relations and provides specific 



predictions that are experimentally accessible. Possible 
extensions to this work include: (i) detailed comparisons 
with experiments and phenomenological models; (ii) en- 
ergy transfer between a tightly bound degree of freedom 
and a weakly bound one [H [T3] ; (in) multiple-step sys- 
tems [13) : (iv) large tilts where inter-band coupling is 
important; and (v) the inclusion of inertial forces |18j . 
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